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We study an extended Bose-Hubbard model with next-nearest-neighbor (NNN) hopping in a 
shaken optical lattice. We show how mean-field phase diagram evolves with the change of NNN 
hopping amplitude t 2 , which can be easily tuned via shaking amplitude. As t 2 increases, a Z 2 - 
symmetry-breaking superfluid (Z 2 SF) phase emerges at the bottom of the Mott lobs. The tricritical 
points between normal superfluid, Z 2 SF, and Mott insulator (MI) phases are identified. We further 
demonstrate the tricritical point can be tuned to the tip of the Mott lobe, in which case a new 
critical behavior has been predicted. Within random-phase approximation, excitation spectra in 
the three phases are obtained, which indicate how the phase transitions occur. 


I. INTRODUCTION 

Ultracold atoms condensed in periodically shaken op¬ 
tical lattices have shown novel properties. Two kinds of 
lattice shaking techniques have been developed. One is 
the off-resonant lattice shaking, in which the shaking fre¬ 
quency is tuned to be very large compared to the band 
gap and width. The hopping parameters and the in¬ 
terparticle interactions can be tuned by lattice shaking, 
which could result in synthetic gauge fields urn, an effec¬ 
tive attractive Fermi-Hubbard model [5] , or topologically 
nontrivial phases ISIEI. 

The other is the near-resonant lattice shaking, in which 
the shaking frequency is tuned to be a little larger than 
the gap of two energy bands. In this case different Bloch 
bands are hybridized, which will dramatically modifies 
the single-particle dispersion and leads to interesting 
phenomena. In a shaken one-dimensional optical lat¬ 
tice, a Z 2 -symmetry-breaking superfluid (Z 2 SF) phase 
has been observed [8], and an effective field theory has 
been constructed to study the normal superfluid-(NSF- 
)Z 2 SF-Mott insulator (MI) phase transition [9]. And the 
effective theory predicted a new critical behavior nearby 
the tricritical point with particle-hole symmetry in three 
dimensions [9]. Algebraical orders nnum and topologi¬ 
cal nontrivial phases [12] are predicted in shaken higher¬ 
dimensional optical lattices. 

The Bose-Hubbard (BH) model, which consists of 
nearest-neighbor (NN) hopping and on-site interaction, is 
used to study a Ml-superfluid transition m- The model 
is a good approximation in the tight-binding limit and 
has been realized in an optical lattice [T4|. Consider¬ 
able efforts have been dedicated to extend the model by 
adding terms, such as next-nearest-neighbor (NNN) hop- 
ping [13 [TS], nearest-neighbor interaction m, dipolar 
interaction m, interaction-induced hopping term m, 
spin structure [20|, or disorder [T3[ [21]. The NN and 
NNN hopping parameters can be renormalized in a dif¬ 
ferent way by off-resonant lattice shaking, hence the ratio 
between them can be tuned [T6| . 

In this paper, we show an extended Bose-Hubbard 
(EBH) model with NNN hopping can be easily realized 
by shaking optical lattices resonantly. Within mean-field 
theory, we find NSF, Z 2 SF and MI phases. We further 


show Z 2 SF phase emerges at the bottom of the Mott 
lobes for nonvanishing NNN hopping amplitude. In three 
dimensions, a new critical exponent of superfluid transi¬ 
tion is predicted near the tricritical point with particle- 
hole symmetry [9]. Nevertheless, the analysis is based 
on a constructed effective theory, and the existence of 
the particle-hole-symmetric tricritical point is in doubt. 
Here within the microscopic EBH model, we demonstrate 
the tricritical point always exists and can be tuned to the 
tip of a Mott lobe. This makes previous work |9] more 
reliable. In the end, we calculate excitation spectra in 
the MI and superfiuid phases in the random-phase ap¬ 
proximation. We find gapless superfluid excitation has 
a quadratic dispersion near the condensate momentum 
at the NSE-Z 2 SE transition boundary. We also demon¬ 
strate in the Z 2 SF phase that the excitation spectrum 
has a roton structure in the strong-coupling limit. 

II. MEAN-FIELD PHASE DIAGRAM 

Let us consider a Chicago-type experiment |H|. Two 
counterpropagating laser beams are time-periodically 
modulated. The Hamiltonian |8] reads 

H{t) = + V COS^{krX + ^^), (1) 

2m 2 

where hkr is the photon momentum, 0{t) = fcos{ujot), 
f and ujQ are shaking amplitude and frequency, respec¬ 
tively, and A = f/(2kr) is the maximum displacement of 
the lattice. 

By performing a transformation, x ^ x — Acos(co’ot), 
in the comoving frame the Hamiltonian reads 

H{t) = ^ + v cos^ {krx) - TlfA ^ (2) 

2m m 

where Ax(t) = mcJoAsin(co’ot). An ac electric field 
Ex = — rnc^o^ cos(ci;o^) is effectively imposed to bosons 
condensed in the unshaken lattice. The first two terms, 
representing the unshaken lattice, give a static band 
structure €x{kx) with Bloch state We choose 

the Bloch state as a basis in our following analysis. 

In the experiment [8], shaking frequency is tuned to 
make s and p bands near-resonant. So we will use two- 
band and rotating-wave approximations in the following 
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analysis. A quasienergy spectrum obtained by numeri¬ 
cally diagonalizing Floquet operator Texp{-~ 
for the lowest 21 bands [8] is shown in Fig. [^(b), where 
T denotes time ordering and T = 27r/co’o is time period. 
Fig. 0 (b) indicates the approximations are very good. 

The Hamiltonian in the tight-binding form reads 


"«) = E (*U.. 

kx 



where 

f —4hp sm{kxd) sin(co’o^) ^ 

I 2iQkx sin(co’o^) — sm{kxd) J ’ 

^kx — ^sp T hgp\ cos(^kxd^^ 

hsp = -^^^{wp{x)\ipx\ws{x)), 
hsp = - ^^^^{wp{x)\ipx\ws{x - d)), 


(4) 

(5) 

(6) 

(7) 


^ is creation operator of a boson in the A-band with 
quasimomentum kx^ d = it jkr is the lattice constant, wx 
is the Wannier function for the A band, A denotes s or p, 
and (... I ... I ...) denotes a real-space integral j dx • • •. 



FIG. 1: Band structure with lattice depth V = 7Er and de¬ 
tuning 5 = 0.44Fr, where Er = {TikxY/{2m) denotes photon 
recoil energy, (a) Band structure before shaking. The ma¬ 
genta and blue lines denote s and p bands, respectively. The 
red line denotes the dressed s band with a detuning 6. (b) 
Quasienergy dispersion of the upper hybridized band. The 
red and green lines are calculated using 2 and 21 bands, re¬ 
spectively. The upper and lower lines with the same color 
denotes the dispersion for shaking amplitude / = 0.2 and 
/ = 0, respectively. 

In the rotating-wave approximation, the effective 
Hamiltonian reads 

^piX'x) 

^kx ^sij^x) T klUJQ 

where 

^ ( 0 ) ■ 




Here we neglect the fast-rotating terms. Before lattice 
shaking, the s band is decoupled with the p band due 
to inversion symmetry (IS). We notice lattice shaking 
effectively breaks IS, causing the coupling between s and 
p bands. Lattice shaking plays the same role with the 
electric field applied in the orbital Rashba effect [22] . 

The quasienergy spectrum calculated by diagonalizing 
the effective Hamiltonian in Eq. <§ is shown in Fig. 
Before shaking, the dressed s band has a perfect cos{kxa)- 
type dispersion, and therefore NNN hopping can be ne¬ 
glected. Lattice shaking changes the dressed s band into 
a hybridized band, in which bosons will stay when turn¬ 
ing on shaking adiabatically. As shaking amplitude in¬ 
creases, the hybridized band dispersion deviates from the 
cos{kxd) form. So an extra cos{2kxd) term needs to be 
considered. 

Assuming the ground state is with quasimo¬ 

mentum kc^ which breaks Z 2 symmetry spontaneously 
for nonvanishing kc [9], time-average interaction energy 
reads 


emt(fcc) = ^ y dtg J dx\^kA^)\‘^^ (9) 

where g is the repulsive interaction strength. 



FIG. 2: Lattice shaking induced NN and NNN hopping. Pa¬ 
rameters {V/Er,S/Er) for dotted, solid, and dashed lines are 
(8,0.25), (7,0.25), and (7,0.45), respectively. Interaction en¬ 
ergy is gn = O.lEr, where n denotes particle density. The red 
and blue lines denote ti/U and t 2 /F, respectively, where U 
denotes on-site interaction energy. 


We project Hilbert space into the upper hybridized 
band. In the tight-binding limit, the next-next-nearest- 
neighbor hopping strength is smaller than the NNN hop¬ 
ping strength and will be neglected without affecting the 
following qualitative results. The off-site interaction en¬ 
ergy is much smaller than the on-site interaction energy 
and will also be neglected. So we will study an EBH 
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Hamiltonian 


Hebh = -ti ^ alaj+t2 ^ a\aj - i 


<«,j> 


«i,j» 


U 


-^ni(ni - 1), 


( 10 ) 


where is the boson operator annihilating boson at site 
rii = a\ai is boson number operator, fi is the chemical 
potential, U is the on-site interaction, and the summa¬ 
tions for the first and second terms are over NN and NNN 
sites, respectively. 

Within standard mean-field theory, the EBH Hamilto¬ 
nian in Eq. (10) can be rewritten as 


HEBH = HMF—ti ^ ataj+t 2 ^ i 
<i,j> «i,j» 

where 


d; 


is needed for a decreased V to 
Eor a fixed E, 


(15) 


V. So an increased 
meet the transition condition in Eq 
as S increases, 2 remains the same, the change of 
with respect to / gets slower, and hence fc increases. 

We numerically calculate the order parameter in 
Eq. (13) by the standard self-consistent approach and 
find the Ml-superfiuid transition is a second-order tran¬ 
sition. So we can use the Landau theory m of phase 
transitions. By using perturbation theory near the MI 
phase boundary for small - 0 , we obtain the mean-field 
ground state energy 

- 1 ) + [1 - 2ix{n, n)] 

(16) 


(11) where M is the number of lattice sites and n) = 


n+l 


Un—fi fi—U(n—l)' 

boundary is given by 


So the Ml-superfluid transition 


( 12 ) 

i 

Hmf = 2 #^ - + i’iO'i) - 

+ - 1), (13) 

hi = ai — ipi represents fluctuation, and t = ti cos{kxcd) — 
t 2 cos{2kxcd). The order parameter pji = (a^) = 
is site-dependent, where (• • •) denotes expectation value 
in the mean-field ground state, kxc is the condensate mo¬ 
mentum in the x direction, Xi is the x coordinate of the 
ith lattice site, and 7/1 is positive and uniform. The mean- 
field Hamiltonian Hme breaks U(l)xZ 2 symmetry when 
pj and kxc are nonvanishing. There are three possible 
phases: (1) MI phase with '0 = 0, (2) NSE phase with 
-0 7 ^ 0 and kxc = 0, and (3) Z 2 SE phase with '0 7 ^ 0 and 

^xc 7^ 0. 

By minimizing the single-particle dispersion e{kx) = 
—2ticos{kxd) + 2 t 2 Cos{ 2 kxd) with respect to kx^ in su¬ 
perfluid phase (0 7 ^ 0 ), kxc has the value of 

r 0 , ti>4|t2i 

kxc=l aarccosj^, \ti\ < A\t2\ . (14) 

I f > < - 4 |i 2 | 

The critical shaking amplitude fc of the NSE-Z 2 SE tran¬ 
sition is determined by the condition 

ti=A\t2\. (15) 


l-2ix{li/U,n)=Q. (17) 


And the boundary condition can be rewritten as 



0 0.05 0.1 0 0.05 0.1 0.15 0.2 



0 0.05 0.1 0.15 0.2 0 0.1 0.2 0.3 0.4 0.5 

ti/U ti/U 


Fig. [2] shows that as shaking amplitude increases, NN 
hopping parameter ti decreases, while NNN hopping pa¬ 
rameter t 2 increases. Here the on-site interaction energy 
U is almost a constant for small shaking amplitude /. 
Eor a fixed detuning, as lattice depth V decreases, initial 
before shaking increases. Initial before shaking is 
almost vanishing. When the detuning 6 is fixed and / 
is small, the degree of the hybridization and the slope of 
the line ti^ 2 -f at any / are nearly the same for different 


FIG. 3: Phase diagram for different t 2 - ^2 increases from (a) 
to (d). The green, red, and blue regions denote the regions of 
MI, NSF, and Z 2 SF phases, respectively. 

The mean-field phase diagram for a fixed t 2 is shown 
in Fig. As ^2 increases from zero, Z 2 SF phase be¬ 
gins to appear at the bottom of Mott lobes near in¬ 
teger values of ja/U. Tricritical points lie on sides 
of Mott lobes. For a fixed filling number n, the tip 
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of the Mott lobe lies at chemical potential {ii/U)c = 
yjn? n — 1, which is the same as that in the stan¬ 
dard BHM. The Z 2 SF region grows and the Mott lobe 
gets thinner and longer because of competition between 
NN and NNN hopping. For a critical NNN hopping 
amplitude (t 2 /t/)c = l/[bb^x((/i/b^)c 5^)]5 the tricritical 
point coincides with the tip of the Mott lobe. When 
t 2 continues to increase, the Mott lobe gets first longer 
and then shorter and finally vanishes. Here the mi¬ 
croscopic theory supports the existence of the tricriti¬ 
cal point in Ref. [9]. The Mott lobe is the longest for 
{h/U)i = l/[4[/x((/i/[/)c,n)], which is larger than 
The critical behavior near a particle-hole-symmetric tri¬ 
critical point is usually different from the mean-field re¬ 
sults and attracts a lot of interest. In three dimensions, 
a 0(2) rotor universality class [11] and a new universal¬ 
ity class |9] have been predicted. For a given n, one can 
tune the parameters to [(/i/l/)c, {ti/U)c, {t 2 /U)c], where 
{h/U)c = 2/[3Ux{{fJ^/U)c^n)]^ to make the tricritical 
point meet the tip of the Mott lobe. 

The Mott lobes have varying shapes and fixed chemi¬ 
cal potentials for their tips for different NNN hopping in 
the mean-field level. A beyond-mean-field theory using 
a 1/(1) quantum rotor approach has predicted the same 
result of bosons with NNN hopping in a two-dimensional 
square lattice ng. The approach only describes the 1/(1) 
symmetry-breaking NSF-MI phase transition. 



Ml(n=1) 


0 0.2 0.4 

Shaking amplitude f 

FIG. 4: Mean-field phase diagram with lattice depth V = 
7Er and shaking frequency loq = bAErlh. l/o is the on-site 
interaction energy before shaking. 


Fig. Il shows the phase diagram expressed in / and 
/r terms. When the shaking amplitude / increases, ti 
decreases and I 2 increases. So an initial NSF phase can 
turn into a Z 2 SF phase, and the tricritical point can be 
tuned onto the tip of a Mott lobe. The parameter regime, 
where tricritical point at the tip of the n = 1 Mott lobe 
lies, is shown in Fig. Figure ^ (a) shows the relation¬ 
ship between critical shaking amplitude fc and lattice 




FIG. 5: Parameters regime for a tricritical point at the tip of 
the n = 1 Mott lobe. 


depth V (detuning (5) as discussed before. We know crit¬ 
ical {ti/U)c is a constant in the parameter regime. From 
previous analysis we also know G at the NSF-Z 2 SF tran¬ 
sition boundary increases as V decreases and does not 
change much for small S. And U is proportional to inter¬ 
action strength g and changes little with V and S. So the 
critical interaction strength Qc in the parameter regime 
increases as V decreases and changes little with 6 un¬ 
der the near-resonant condition, as shown in Fig. |5](b). 
When a shaken one-dimensional lattice system is tuned 
to this tricritical point, one can measure the critical expo¬ 
nent via the in situ technique [24| and a new universality 
class is expected |9]. 


III. COLLECTIVE EXCITATIONS 

In this section, we study collective excitations at zero 
temperature. Following the standard-basis operator ap¬ 
proach [25j[26], we choose eigenstates ||i(a) | of the single¬ 
site mean-field Hamiltonian in Eq. ( [l^ as a basis, 

and the EBH Hamiltonian Hebh in Eq! ^11[ ) can be 
rewritten as 

Hebh = T + (-^i E +^2 N ) 

i,cx <hj> <<i,j>> 

aa' gp' 

where E^ is the mean-field energy per site, = 

= {ia\al\ia'}{jl3\aj\jl3'). 

The single-particle retarded Green’s function is defined 
as 

= -ie{t - t'){[ai{t),al{t')]), (20) 

where 0(t) is the step function. In the standard basis, 
the Green’s function reads 

9i,j{t — t)= ^ “ O 5 ( 21 ) 

aa' gg' 

where 
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By introducing the random-phase approximation, one 
obtains the equations of motion for G in the frequency 
and momentum space, 

^aP'^a'P^aa' — )^aa'/3/3 '5 

Daa' ^ ^ “ 1 “ ^xc)Tcx'a,^^' 

77 ' 

“1“ a'a ^77'/3/3' {^x; ^); 

(23) 

where Daa' = (Laa) - {La'a'), Kc IS given in Eq. 


(14) for superfluid phase and is zero for MI phase, 
Ta'a^y = Va'a = (fa'|e*''-="*a||zQ;), and = 

Ta^aj'y' is site independent. Equa¬ 


tions (23) are linear equations of Uaa'Gaa'(3^' and 


Z)c««' yIl'^aa'/ 3 / 3 '- Substituting the solution into the 
Green’s function g{kx,uj)^ one obtains 


/7 \ _ ^{kx 2kxc’>^) 

“ l-e{kxMkx-2kxc.ujy 


where 


A yOaVio VaOVoa 

,;+-AE„“a;+ + AE„ 

a L 

^12 (<^) = ^21 (‘^) 

= E 


(24) 

(25) 

(26) 



VaOVOa 


(27) 

— AEa 

uj^ + AE(y _ 

") 

VOaVaO 

yoL^y^oc 


(28) 

— AEc^ 

uj-\- + AEq, 

5 

2 mean-field single-site ground 

state. 


^22 (^) - ^ 


co’+ = cj + i0+, and AE^ = — Eq. Here the Green’s 

function is a generalization of that in the standard BHM 
[26], in which bosons condense at zero momentum {kxc = 
0 ). 

In the MI phase, the basis is just the Fock state \ia = 
n). For a commensurate filling n, there is nonvanishing 
Vn'n = V^^n',n-i- The Green’s function reads 


g{kx,uj) = 


1-Z 


where 


^p,h — 


CJ+ - Ep uj^- Eh' 
U {2n — 1) — 2/i + €{kx) 


(29) 


Z = 


± + 2(2n + l)Ue{kx) + e(A)^)2 

g + U E Ep(kx) 

Ep{kx) Eh{kx) 


(30) 

(31) 


Ep^h represents particle (hole) excitation. The condition 
of existence of a gapless excitation at kx = kxc exactly 
gives the MI phase boundary in Eq. (18). 




Quasimomentum Quasimomentum k^kr 

FIG. 6: Excitation spectra with parameters marked in Fig. 
(a) and (b). Red, blue (dashed), and black lines denote ex¬ 
citation spectra in the superfluid phase, at the Ml-superfluid 
transition boundary, and in the MI phase, respectively. 


In the superfluid phase, we will numerically calculate 
the Green’s function in Eq. ( [^ and the spectral func¬ 
tion A{kx^03) = — (I/tt) Im g(^,uj), of which excitation 
modes give excitation spectra. Fig. shows excitation 
spectra near different Ml-superfluid phase boundaries. 
There are two gapless spectra in the superfluid phase with 
positive and negative energy corresponding to quasipar¬ 
ticle and quasihole excitation, respectively. In the super¬ 
fluid there are also gapped excitation modes as a conse¬ 
quence of the band structure of the lattice system. In the 
Z 2 SF phase, roton excitation spectrum has been observed 
in the weakly-interacting regime m- Here we show the 
roton excitation spectrum in the strong-coupling limit 
in Fig. [^(a). Figures]^ (a) and (c) show linear disper¬ 
sion around condensate momentum in the Z 2 SF and NSF 
phase, respectively. Figures]^ (b) and (d) show quadratic 
dispersions around kx = 0 dit the NSF-Z 2 SF transition 
boundary without and with particle-hole symmetry, re¬ 
spectively. In the superfluid phase, the quadratic dis¬ 
persion indicates stronger phase fluctuations and weaker 
superfluidity than linear dispersion HU. At the lower Ml- 
superfluid transition boundary, the gapless particle dis¬ 
persion vanishes, the hole dispersion becomes quadratic 
around kx = kxc^ and a Mott gap is opened, which means 
disappearance of superfluidity. At the lobe tip, the Mott 
gap vanishes. In the Mott phase, the dispersions of both 
particle and hole excitations are gapped. 
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IV. CONCLUSION 

In conclusion, we have shown a significant NNN hop¬ 
ping effect in near-resonantly shaken optical lattices. 
We studied the mean-field phase diagram for a one¬ 
dimensional EBH model and found tricritical points be¬ 
tween three phases. Furthermore, we calculated corre¬ 
sponding microscopic parameters to the EBH model pa¬ 
rameters and provided strong support for the existence 
of the tricritical point with particle-hole symmetry. A 
new critical behavior [9] is expected to be verified by the 
in situ technique [24] in the parameters regimes. We 


also calculated the excitation spectra in all three phases 
and showed how the spectrum evolves during the phase 
transitions. In the Z 2 SF phase, the excitation spectrum 
has the roton structure. At the NSF-Z 2 SF transition 
boundary, the quasiparticle (quasihole) excitation has a 
quadratic dispersion relation around = 
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